arXiv:astro-ph/0112434vl 18 Dec 2001 


Draft version February 1, 2008 

Preprint typeset using IAT^X style emulateapj v. 14/09/00 


A NEW WAY TO DETECT MASSIVE BLACK HOLES IN GALAXIES: THE STELLAR 

REMNANTS OF TIDAL DISRUPTION 

R. Di Stefano 

Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138 
Department of Physics and Astronomy, Tufts University, Medford, MA 02155 

J. Greiner 

Astrophysical Institute Potsdam, 14482 Potsdam, Germany 
S. Murray, M. Garcia 

Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138 
Draft version February 1, 2008 

ABSTRACT 

We point out that the tidal disruption of a giant may leave a luminous (10 35 —10 39 erg/s), hot (10 — 100 
eV) stellar core. The “supersoft” source (SSS) detected by Chandra at the center of M31 may be such a 
core; whether or not it is, the observations have shown that such a core is detectable, even in the center 
of a galaxy. We therefore explore the range of expected observational signatures, and how they may be 
used to (1) test the hypothesis that the M31 source is a remnant of tidal stripping, and (2) discover 
evidence of black holes and disruption events in other galaxies. 

Subject headings: X-rays: stars — galaxies: M31 


1. TIDAL DISRUPTION OF GIANTS: THE CORE REMAINS 


A star of mass M* and radius i?* can be tidally dis¬ 
rupted by approaching within R t of a black hole (BH). 


Rt = 



(1) 


Myh is the BH’s mass; 77 is a parameter of order of unity. 
Let R s be the BH’s Schwarzschild radius. For M^h ~ 
10 s Mq, R t < Rs for most main-sequence stars, and only 
tidal disruptions (TDs) of giants lead to observable effects. 
The rate of such TDs, Nq, depends on Mbh, the kine¬ 
matics of stars near the galactic nucleus, and the mass 
and age distribution of these stars (see, e.g., Magorrian 
& Tremaine 1999). Computed rates are typically in the 
range 1CD 6 — 10~ 4 yr _1 . A convenient expression is 
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with L* = 1.8 x 10 lo Z/ Q ; this was derived by Syer & Ulmer 
(1999), who also provide cautions about the use of a sim¬ 
ple scaling law. Studies of TD have focused on possible 
associated accretion events, lasting for months or decades, 
with luminosities as high as ~ 10 44 — 10 46 erg/s (see, e.g., 
Hills 1975, Lidskii & Ozernoi 1979, Gurzadyan & Ozernoi 
1980, Rees 1988, Loeb & Ulmer 1997, Ulmer, Paczynski, 
& Goodman 1998, Ulmer 1999). 

The Hot Core of a Disrupted Giant: The disruption 
leaves an end-product, the giant’s hot dense core, whose 
presence, influence, and observability are the subjects of 
this paper. The core remains hot (T > 10 5 K) and bright 
(L > 10 2 Lq) for 10 3 — 10 6 years, thereby providing the 
longest-lasting signal of a TD. 

Detectability: The soft X-ray sensitivities of Einstein, 
ROSAT, Chandra, and XMM-lVewton have allowed them 
to detect supersoft X-ray sources (SSSs) in other galaxies. 
SSSs (Greiner etal. 1991, Rappaport, Di Stefano & Smith 
1994, Kahabka & van den Heuvel 1997) have luminosities 


(10 35 — 10 38 erg/s) and temperatures (k T ~ 10 — 100 eV), 
comparable to those of hot stellar cores. Chandra’s 0.5-1" 
angular resolution allows individual sources to be resolved, 
even in the dense central regions of nearby galaxies. 

The Example of M31: We have studied 8.8 ksec of data 
collected by Chandra’s ACIS-I detector. To identify SSSs 
we applied a hardness ratio test to all point sources in a 
16' x 16' region centered on the center of M31, selecting 
all sources with more than 50% of their photons below 0.7 
keV. Three sources satisfied this condition. Remarkably, 
one of these is coincident with the center of M31. The lu¬ 
minosity is ~ 10 37 erg/s; the data are consistent with little 
or no emission above 1.5 keV. The source appears to be 
variable; data from earlier missions ( Einstein , ROSAT), 
e.g., are consistent with the source providing zero flux on 
one occasion. Garcia et al. (2000) have argued that the 
density of point X-ray sources in this field is so low that a 
chance coincidence with the center is highly unlikely; this 
argument is even stronger when applied to SSSs. The ob¬ 
served SSS-like behavior is therefore most likely related to 
the environment of the center, presumably to the presence 
of a MBH. The data do not seem consistent with ADAF 
disk models (see, e.g., Garcia et al. 2000), but are com¬ 
patible with the signature of a hot stellar core. Whether 
or not the source is a hot core, the M31 observations es¬ 
tablish that it is possible to detect such an object in the 
centers of nearby galaxies. 

2. THE NATURE OF THE CORES 

The nature and characteristics of the core depend on 
the initial mass, M* of the stripped star, and on its state 
of evolution at the time of disruption. 

2.1. M, < 8 M q 

Tidal stripping brings a premature halt to nuclear burn¬ 
ing. While any residual nuclear burning continues, the lu¬ 
minosity remains at the value it had when the star was a 
giant. The radius shrinks, and the effective temperature 
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increases. Depending on the core mass, these systems will 
stay at L > 10 2 Lq, and T > 10 5 K for 10 3 — 10 5 yrs (see, 
e.g., Blocker 1995; Gorny, Tylenda & Szczerba 1994). 

Nature has provided direct analogs, in the form of the 
hot cores at the centers of planetary nebulae. A notable 
example is IE 0056.8-7154, a planetary nebula system in 
the Small Magellanic Cloud, with L = 2 x 10 37 erg/s and 
k T = 30 - 40 eV (Wang & Wu 1992). 

The evolution of the hot core is influenced by its mass, 
its state at the time of disruption (e.g., its thermal pulse- 
cycle phase), and the mass of the original star (see, e.g., 
Iben 1995 and references therein). The hot core of a star 
that was disrupted at some arbitrary time in its evolu¬ 
tion, is descended from a more massive star than a core 
of the same mass emerging during the final stages of an 
uninterrupted evolution. The core of the disrupted star is 
therefore younger, and is generally less compact, with a 
significantly higher core temperature; it cools more slowly 
than typical post-AGB stars. (See, e.g., Blocker 1995; 
Gorny, Tylenda, & Szczerba 1994.) Although these re¬ 
sults were arrived at by studying winds, they are expected 
to be valid for the extreme mass-loss scenario associated 
with disruption. 

The stellar mass function favors stars in this mass range, 
and they are therefore expected to be disrupted more fre¬ 
quently than stars of higher mass. Nevertheless, we go on 
to consider high-mass stars, because their disruption may 
also produce remnants with observable signatures. 

2.2. M* > 8 Mq 

The stripped core is a helium star (He star). Two pro¬ 
cesses (in addition to TD) may produce He stars. First, 
a very massive star may eject such copious winds that it 
effectively ejects its own hydrogen envelope. Second, a 
massive star in a binary may overfill its Roche lobe as 
it evolves and lose its hydrogen envelope by either trans¬ 
ferring it to its companion (if the mass ratio is not too 
extreme), or losing it during a common envelope episode. 
The evolution of isolated He stars and He stars in binaries 
has been well-studied; see, e.g., Maeder (1981), Delgado 
& Thomas (1981), Habets (1986, 1987), de Loore & De 
Greve (1992), Reese (1993), Woosley, Langer, & Weaver 
(1995). The characteristics of each stage of the He star’s 
evolution depend on many factors, including the mass of 
the star when its envelope was lost, and the continuing 
rate of mass loss to winds. A common feature of these 
evolutions is, however, that for a wide range of system 
parameters, there is a time interval lasting up to several 
million years, during which the star is both very luminous 
(10 38 — 10 39 ergs/s) and very hot, with kT near 10 eV. 
Characteristic evolutionary sequences include an epoch of 
helium burning, a shorter epoch during which heavier ele¬ 
ments are burned, and, in some cases, a Type lb or Type 
Ic, hydrogen-poor, supernova. 

Typical temperatures are lower, and typical luminosities 
higher than those considered in §2.1. Another distinction 
is that the core may not be hottest and brightest immedi¬ 
ately after the TD. This is because the further evolution 
of the He star will typically increase its luminosity and 
alter its temperature; in addition, winds that can block 
emergent soft X-rays tend to decrease with time. 

Wolf-Rayet stars may be the closest well-observed 


analogs of the He-star remnants of TDs. They tend to eject 
such significant winds, however, that their appearance may 
be dominated by the interaction of this surrounding mat¬ 
ter with radiation from the He star. Perhaps because of 
this, no known He star has a “supersoft” signature. It is 
possible, however, that some of the known SSSs, whose 
fundamental nature is not yet understood, could be He 
stars. 

He-star remnants of TDs may be more likely than other 
He stars to be observed as SSSs. This is because the 
BH continues to act as a sort of vacuum cleaner for mass 
outside the tidal radius (~ R t ,]$). Furthermore, 

since self-generated winds tend to decrease with time (see, 
e.g., Woosley, Langer, & Weaver 1995), there will typically 
be less and less material within the tidal radius, even as 
nuclear evolution continues. Thus, He stars near MBHs, 
especially those that have undergone TD, may provide the 
cleanest signatures of the underlying properties of the star, 
distinct from the signatures of its interactions with its en¬ 
vironment. 

3. OBSERVATIONAL SIGNATURES 

Luminosity: For both white dwarfs and He stars, the 
luminosity of the stripped core should not exceed that of 
typical giant stars. Luminosities in excess of ~ 10 39 erg/s 
must be due to some alternative or additional effect. 
Spectrum: The maximum temperature expected for hot 
white dwarfs (He stars) is ~ 100 (15) eV. Therefore, when 
Nh is low, the X-ray spectrum is dominated by soft pho¬ 
tons. There can also be a modest contribution from pho¬ 
tons at higher energies. First, the spectrum could differ 
from a simple thermal shape, e.g. due to opacity effects, 
as is well known in the case of white dwarf atmospheres 
(see, e.g., Dupree & Raymond 1983, van Teeseling et al. 
1994, Hartmann & Heise 1997); thus, atmospheric effects 
must be included when fitting the X-ray spectrum. Sec¬ 
ond, there could be contributions from a corona, an infall 
disk within the stellar core’s Roche lobe, or shock heating 
of the ambient medium. Finally, it may not be possible 
to resolve the hot core and the MBH, which should be 
accreting matter at a low level, e.g., from stellar winds. 
Note that any radio emission is more likely to be originat¬ 
ing from the region around the MBH than from the hot 
remnant of TD. 

It is important to look for distinctions between the signa¬ 
ture of tidally stripped remnants and other possible ex¬ 
planations for SSS signatures. One distinction, e.g., be¬ 
tween SSSs in binaries and stripped cores, should be that 
the distinctive disk signatures of the former (Popham & 
Di Stefano 1996) should be missing in the latter. 
Ionization: Even if the total flux from a stripped core 
over its lifetime (which, for the purposes of ionization sig¬ 
natures, may be defined to be the time interval during 
which T > 50, 000 K) is smaller than that from the accre¬ 
tion flare event, the number of high-energy ionizing pho¬ 
tons it emits may be significantly larger. The distinctive 
signature is photoionization of high-energy states. 

Time Variability: Intrinsic variability may occur at 
many time scales. Here we concentrate on variability due 
to absorption of soft X-radiation by a varying amount of 
matter lying between the source and our detectors. Figure 
1 demonstrates that the decline in count rate with increas- 
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ing column density can be dramatic. 



Figure 1: Chandra ACIS-I detector count rates of a black 
body radiator with a bolometric luminosity of lCP 8 erg/s in 
M 31 for different temperatures (solid lines) and depending 
on the intervening absorbing column. 


If the stripped core orbits the MBH, the flux we receive 
may decrease as it travels behind gas in the disk around 
the MBH. If the mass and radius of the star were M* and 
R * before the TD, and if the orbital separation between 
the core and the MBH is presently equal to the separation 
at which the star filled its Roche lobe, then 


Porb = 0.32 years ry 


R * 
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If the core escapes from the vicinity of the MBH, it 
could also exhibit variability as it passes behind gas clouds. 
Making a set of simple assumptions, we estimate the 
amount of time, r, it takes for the hot core to pass be¬ 
hind a cloud with column density Nh, located a distance 
a from the massive black hole . 1 
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where £ is of order unity, U 150 is the transverse speed in 
units of 150 km/s, and Nh,o = 1.9 x 10 22 . The actual 
duration of obscuration may be shorter or longer, depend¬ 
ing on the geometry of the cloud and the track of the hot 
remnant. 

Doppler shift: The orbital speed of the star in a cir¬ 
cular orbit around a 3 x 1O 7 M 0 black hole, with orbital 
separation, a, such that a 100 Rq star fills it Roche lobe, 
is ~ 5% c; the speed of the hot remnant could be greater. 
Thus, the most direct confirmation that an observed SSS 
is a tidally stripped remnant, would be a periodic Doppler 
shift compatible with an orbit small enough to allow TD 
of the progenitor star. Significant blue or red shifts could 
possibly be inferred from spectral fits. Even if the remnant 
does not orbit the MBH, it might be possible to measure 
its velocity and establish the likelihood that it is escaping 
from the central gravitational potential. 


Spatial Distribution: If the hot stellar core escapes 
from the region around the MBH , traveling along a direc¬ 
tion perpendicular to our line of sight, it could traverse a 
considerable distance before it dims and cools. If the av¬ 
erage speed is ~ 10 8 cm/s, the core could travel ~ 1 ( 0 . 1 ) 
kpc in 10 6 (10 5 ) years. If the product of the TD rate and 
average hot-core lifetimes is larger than unity, we might 
expect to see several hot cores ringing the MBH, with the 
oldest cores further away. Any cores located ~ 1 kpc from 
the galaxy center would likely be He stars. It is interesting 
to note that the other two SSS candidates near the center 
of M31 are each several arcminutes from the galaxy’s cen¬ 
ter, produce count rates similar to that of the central SSS, 
and are each considerably softer than the central source. 
Population Signatures: (1) Consistency with the 

model: If N M31’s SSSs are remnants of TD events, 
each represents a much larger population of giants in the 
galaxy’s central region at the time of the TD. We must 
therefore test whether the present-day population could 
have evolved from one with the requisite number of giants 
some 10 3 — 10 6 years ago. If we determine that a subset of 
the TD remnants are He stars, it is likely that there were a 
large enough number of centrally-located high-mass stars 
active in the recent past to lead to some supernovae. We 
can ask if the characteristics of the central region (up to 
~ 1 kpc from the center) are consistent with the required 
number of Type II or Type lb and Ic supernovae. The de¬ 
gree and characteristics of the ionization of the central re¬ 
gion must also be consistent with the effects expected from 
the N TD remnants and the requisite total population of 
high-mass stars. It may be instructive to consider the pop¬ 
ulation of He stars near the center of the Milky Way, which 
is much better studied. This population is large and in¬ 
cludes members with apparently unusual properties. (See, 
e.g., Tamblyn et al. 1996, Najarro et al. 1994 and ref¬ 
erences therein.) Unfortunately the large column density 
toward the center of the Milky Way makes it impossible 
to study any soft X-ray emission from these stars. 

(2) Eliminating other possibilities: Population studies 
of the region around the galactic center can lead to esti¬ 
mates of the numbers of post-AGB stars and also of X-ray 
binaries of various types. The physical processes that give 
rise to each of these objects also may give rise to SSSs 
that have nothing to do with TD. We therefore need: (1) 
observations of the central population, ( 2 ) theoretical esti¬ 
mates of the number of SSSs not due to TD expected, and 
(3) comparisons with studies of populations located away 
from the central region. These should allow us to estimate 
the probability that soft sources observed near the center 
have their genesis through processes that are ordinary and 
expected parts of stellar and binary evolution. 


4. OTHER GALAXIES 

Radiation from a hot stripped core can be detected in 
other galaxies if ( 1 ) the count rate is high enough, and ( 2 ) 
the source can be resolved. 

Count Rate: The count rate from the soft X-ray source 
in the center of M31 (Garcia et al. 2000) was ~ 10 


1 Consider a cloud with mass m located a distance a from the MBH. If the cloud is corotating, the largest volume it could have is roughly 
equal to the volume of a sphere whose radius is the effective Roche lobe radius. This allows us to relate the mass of the cloud to its radius R, 
its distance a from the MBH, and to M^h- Assume the cloud is composed of hydrogen and that we can approximate an average value of Nh 
by the number of hydrogen atoms divided by the square of the cloud radius. This allows us to relate R to Nh- Finally, we use r = R/v. 
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counts/ksec. Based on spectral fits to nearby (brighter) 
sources, the absorption was estimated to be 2.8±1.0 x 10 21 
cm -2 . The number of counts expected from an E ksec 
exposure of a source with similar spectrum, located a dis¬ 
tance D from us, is approximately 

/ \ 2 / \ 


Nj = 10 counts x E x 


Mpc 

D 


x/ 


N h 


2.8 x 10 21 


• ( 5 ) 


Nh is the column density between us and the new source. 
For face-on or early-type galaxies Nh could be smaller 
than 2.8 x 10 21 , increasing the count rate in a way that 
depends on the spectrum, but which is likely to be de¬ 
scribed by a non-linear monotonically decreasing function, 
/( 2.8x10^ )• ( See Figure 1.) 

A 50 ksec exposure could therefore yield 20 counts from 
a hot stellar core an in M31-like galaxy as far away as 5 
Mpc; for galaxies with less internal absorption, the count 
rate would be higher. Should the photons’ energies be clus¬ 
tered below 1 keV, we would have convincing evidence for 
a centrally located soft X-ray source. XMM-Newton, with 
its higher effective area, could achieve similar numbers of 
counts for galaxies up to 10 Mpc away. There are about 
200 (600) galaxies within 5 (10) Mpc listed in Simbad and 
the Principal Galaxies Catalog (PGC). The galaxy type or 
orientation could allow the detection of soft central X-ray 
sources in roughly 1/3 of these galaxies. 

Spatial Resolution In M31 there are 4 other X-ray 
sources within ~ 10" of the nucleus. If the central den¬ 
sity of X-ray sources in other galaxies is similar, we might 
not be able to resolve the central source in galaxies farther 
from us than 2 — 3 Mpc. Note, however, that, if the soft 
central source is variable, as the source in M31 appears 
to be, it may be possible to successfully apply to X-ray 
searches for central sources, the same difference imaging 
techniques that have been applied to optical searches for 
supernovae and microlensing. If the hard X-radiation is 
steady, or varies in a way well-known from studying other 
X-ray stars, then we may be able to infer the existence of 


the soft central source. XMM-Newton, which has 4 — 5 
times Chandra’s effective area at low energies, will make 
such differencing studies possible for M31-type galaxies at 
least as far away as 5 Mpc. This is particularly so if si¬ 
multaneous observations with Chandra can resolve some 
sources and help to determine their spectral characteris¬ 
tics. 

Statistical Tests: Statistical questions can be posed 

when we study a population of galaxies. For example, are 
galaxies for which there is good evidence (perhaps from dy¬ 
namical studies) for a MBH more likely to house a centrally 
located population of soft X-ray sources, than galaxies for 
which the dynamical studies do not indicate the presence 
of a MBH? 


5. CONCLUSION 

When a giant is tidally disrupted, the stripped core may 
provide a hot and bright signature for 10 3 — 10 6 years. 
Detecting such cores and relating them to TD events is 
challenging. Fortunately, the present generation of in¬ 
struments, including Chandra, XMM-Newton, and HST, 
promise opportunities for progress. 

This new method of studying MBHs complements 
searches for flare events that may be due to TDs. The 
latter must of course be used for distant galaxies, and 
there is a growing body of evidence for luminous UV and 
X-ray events lasting for months to years, which may be 
consistent with the accretion of the envelopes of tidally- 
disrupted stars. 2 Such events are so rare, however, that 
we would be lucky to detect one in the near future among 
the small number of galaxies within 5 or 10 Mpc. In nearby 
galaxies, however, the long-lived signatures of the stripped 
cores may allow us to infer that TDs have occurred. The 
two complementary modes of study should help us to bet¬ 
ter understand the frequency and consequences of tidal 
disruption events and to use them as probes of massive 
black holes in the center of galaxies. 
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